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Abstract 

The Lee-Suzuki iteration method is used to include the folded diagrams in 

(2) 

the calculation of the two-body effective interaction v ^ between two nucleons 
in a no-core model space. This effective interaction still depends upon the 
choice of single-particle basis utilized in the shell-model calculation. Using a 
harmonic-oscillator single-particle basis and the Reid-soft-core NN potential, 
we find that overbinds 4 He in 0, 2, and AhQ model spaces. As the size of 
the model space increases, the amount of overbinding decreases significantly. 
This problem of overbinding in small model spaces is due to neglecting effective 
three- and four-body forces. Contributions of effective many-body forces are 
suppressed by using the Brueckner-Hartree-Fock single-particle Hamiltonian. 
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1 Introduction 



Previous calculations of the shell-model effective interaction involved a number of 
uncertainties. The major ones included the choice of the single-particle (s.p.) basis, 
the choice of the starting energy uo in the Brueckner G-matrix [|TJ calculation, the 
neglected contribution from higher-order core-polarization diagrams, and the effects 
of both the real and effective three- and higher-body forces. It was proposed in Ref . 
that the core-polarization diagrams could be eliminated by adopting a no-core model 
space, in which all the nucleons in a nucleus are treated as active. In Refs. |3|, |4j] , the 
no-core approach has been used and satisfactory results are obtained for light nuclei. 

The results of this no-core approach depend on the choice of the starting energy 
for the G-matrix. One may argue that this uncertainty in the choice of the starting 
energy u can be avoided by evaluating an energy-independent interaction employing 
the iteration methods proposed by Lee and Suzuki 0] and Krenciglowa and Kuo [^J. 
Applying this technique, one can sum the folded diagrams to all orders to obtain a 
starting-energy-independent effective interaction ir ff ; . Therefore, with the use of a 
no-core space and including the folded diagrams, the effective interaction obtained is 
subject to only two of the major uncertainties mentioned above, i.e., the choice of the 
s.p. basis and the effect of the neglected effective three- and higher-body forces. 

These two remaining uncertainties are related. With an optimally chosen s.p. ba- 
sis, the contribution from the effective many-body forces could be minimized [[F|. 
Furthermore, both of the uncertainties are related to the size of the model space and 
are expected to diminish as the size increases. Indeed, the effective interaction is 
only introduced with the truncation of the infinite Hilbert space to a finite-size model 
space. Furthermore, one should keep in mind that so-called Q-box diagrams, like the 
one displayed in Fig. la, are not the only source of effective many-body forces. Even 
if the Q-box is restricted to two-body terms, which means the G-matrix for no-core 
calculations, the inclusion of folded diagrams yields effective three-body forces (see 
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Fig. lb). It has been demonstrated that such many-body forces are non-negligible 
||, in particular if a large number of active particles has to be considered 0. As we 
will discuss below, one may try to minimize the effects of such many-body forces by 
introducing an appropriate auxiliary field. One may view the present effort as an ex- 
tension of the work of Ref. to the case of realistic NN potentials which capitalizes 
on the results presented in Refs. 0, || . 

In this work, we will use the Lee-Suzuki method [0 to calculate the starting- 

(2) 

energy-independent two-body effective interactions v g ff for no-core, harmonic-oscillator 
(HO) model spaces and study the dependence of the shell-model results obtained for 
4 He with t>gg on the HO basis parameter HQ and the size of the model space. It has 

(2) 



been noticed in a previous work |10| that v g ff tends to overbind light nuclei. Here we 



will show that the overbinding is quite significant when the model space is relatively 
small. We will show that the overbinding problem can be cured by introducing an 
auxiliary field, such as the Brueckner-Hartree-Fock (BHF) approximation. 

After this introduction we will present some details on the evaluation of the energy- 
independent effective two-body force for no-core shell-model calculations in section 
2. Numerical results for the binding energy of 4 He will be presented in section 3. In 
section 4 we will discuss the influence of an auxiliary potential, and section 5 contains 
the conclusions of the present investigation. 



2 Calculation of 

The Brueckner G matrix is calculated according to the following equation: 

G{u) = v 12 + v 12 ® -v 12 , (1) 

u — (ht + h 2 + v 12 ) 

where v± 2 is the NN force for which we will use the Reid-soft-core (RSC) potential 



11], us is the starting energy, Q is the Pauli operator which excludes the scattering 



into the two-particle states inside the model space. For a full no-core NhQ space, we 
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define Q as 

Q = for m + n 2 < N 

= 1 for n 1 + n 2 > N. (2) 

where n, = 2n r (i) + l(i) = 0,1,..., are the principal quantum numbers of the 
s.p. states occupied by the two intermediate-state nucleons in the multiple scattering 
process. 

For the first part of our discussion (see section 4 for an alternative choice) the 
s.p. Hamiltonian in Eq.(jl]) is taken as 

h l = U + Ul = ^ + (~m^ 2 rf - Vo) = U + (uf - Vo) = hf° - V , (3) 

where we use a s.p. potential (uj) that is of the shape of a harmonic oscillator but 
is shifted downward by an amount Vo to make it more realistic. The quantity Vq 
represents the depth of the mean field of the nuclear medium. It is convenient to 
define a shifted starting energy as uj' = (u + 2Vo) and rewrite Eq.(|I|) as 

G{LJ,) = Vl2 + Vl2 u>-( h ™ + hr + v 12 ) Vl2 > (4) 

where h uo = ti + uf° is now a pure HO Hamiltonian. Since we will use the Lee- 
Suzuki iteration method H to take into account the folded diagrams, the resulting 
effective interaction will be independent of the starting energy u' as well as the shift 
Vo. Therefore, no specific choice for the value of Vq needs to be made. 

However, in the case when the folded diagrams were ignored (which is a common 
practice in effective-interaction calculations), one would have to choose a reasonable 
starting energy to minimize the contribution from the folded diagrams. It should then 
be noted that the starting energy u' used in Eq.(^), unlike u, does not correspond 
to the energy E 2 of the initial two-particle state in the ladder diagrams. Rather, it 
is related to (E 2 + 2Vo). When the two nucleons in the initial state occupy bound 
s.p. states, E 2 is negative. But u' ~ (E 2 + 2V ) could very well be positive. In fact, 
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it has been found [|T^, that for the two valence neutrons in 18 O, a value of about 
70 MeV for uj' yielded reasonable G-matrix elements. 

In order to obtain the starting-energy-independent two-body effective interaction 
v e Q, we calculate G(uj') of Eq.(|J) for 11 values of &' ranging from about — 5M1 to 
about 5MX These 11 sets of G matrices are then used to numerically calculate the 
derivatives of G(uj') with respect to uj' to the 9th order. Once the derivatives of G(oj') 

(2) 

are obtained, we proceed with the Lee-Suzuki method to obtain v eS ! . Here we point 
out that the number of iterations needed for convergence strongly depends on the 
value of uj' at which the derivatives are evaluated. It generally exceeds the number 
of derivatives retained in the iteration procedure. 

In Fig. 2, we show the values of the diagonal two-body matrix elements (TBMEs) 
of G(oj') and in the states |(0si/ 2 ) 2 ) j,t with J=0, T=l and J—l, T=0 for a wide 
range of uj'. For these matrix elements we use M7=16 MeV and the Pauli operator is 
defined in Eq.(|^) with N=4. It can be seen from the figure that the matrix elements 



of G{u') decrease (i.e., become more attractive) with increasing to 1 ', while those of v { 



(2) 
cff 



are independent of uj' . 

We also note in Fig. 2 that there is a particular value of uj' for which the matrix 

(2) 

elements of G{uj') are about equal to those of v g ff . This observation is the basis of an 



approximation scheme presented and tested in Ref. JTO] and then used in Ref. 



3 Shell-Model Results 

We perform the matrix diagonalization for the shell-model Hamiltonian 

/ a \ A 3 

H S M = (2** _ Tc - m - + E^cfffX?) + ^Coulomb + K H c.m. ~ ^^)- 
\i=l / i<j ^ 



(5) 



,2 



In the above equation the t{ = Yh/ (2m) are the one-body kinetic energies, T c m = 
Pi) 2 /(2mA) is the cm. kinetic energy and Vbouiomb is the Coulomb interaction. 
The proton and neutron masses are taken to be the same. The last term (with A=10) 
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in the above equation forces the cm. motion of the low-lying states in the calculated 
spectrum to be in its lowest HO configuration. 

In Fig. 3, we plot the calculated ground-state (g.s.) energy of 4 He as a function 
of the HO basis parameter HQ for three model spaces of different sizes, OhQ ( "iV=0" 
curve), 2hQ ("iV=2" curve), and 4M7 ("iV=4" curve). In the OhQ model space which 
consists of only the 0si/ 2 major shell, the g.s. energy begins at -33.9 MeV for hQ=10 
MeV, decreases to a minimum of -43.3 MeV at hQ=22 MeV and then increases to 
-41.6 MeV for M7=28 MeV. These results significantly overbind the g.s. of 4 He. They 
are 5.6 to 15.0 MeV lower than the experimental g.s. energy of -28.3 MeV and are 
9.3 to 18.7 MeV lower than the value of -24.6 MeV obtained in the (nearly exact) 
Green's function Monte Carlo (GFMC) approach fl4| using the RSC potential. 

In this one-major-shell model space, there is a simple way to obtain the above 
results. The g.s. energy of 4 He can be expressed in terms of the effective-interaction 
TBMEs as (Oa = 0s 1/2 ) 

' 3 — ^ ~ r '~,2|„,(2)|n„2\ . /n 2. ,(2) 



= 3 ( -M> 



(0 S >y|0 S 2 )j=o,T=i + (Os>^|0^) J=1 ,T=o , (6) 



where the first term is the kinetic energy (with the cm. contribution subtracted) 
and the second term is the effective-interaction energy. According to Ref . ||1 5|| , we 
know that for this one- dimensional model space, the effective-interaction TBMEs are 
related to the eigenenergies of the Schrodinger equation: 

(/ii + h 2 + v 12 )(pj : T = Ej,t 0j,t (7) 

through 

-2L,(2)|n n 2\ ■■ ~ /3, 



(0s 2 \v ( e s>\0s 2 )j, T = Ej, T - 2 ^-m) , (8) 

where ~hfl is the eigenenergy of the s.p. Hamiltonians hi and h 2 for the 0si/ 2 state. 
For htt=22 MeV, we obtained 



E i = 55. 183 MeV and E x o = 45.892 MeV, 
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which lead to 



and 



Therefore, 



(0s 2 \v^\0s 2 ) 0A = -10.817 MeV 



(0s 2 |u2 ) |0s 2 ) 1 , = -20.108 MeV. 



E gs = 3 Q x 22^ + 3(-10.817 - 20.108) = -43.275 MeV, 

which agrees with the result that we obtained through the Lee-Suzuki iteration pro- 
cedure. We further remark that the effective-interaction TBMEs (0s 2 \v^\0s 2 )j i t are 



equal to Gj^ioo') with 00' = Ej^- This is another property observed in Ref.[0 for 
the effective interaction in a one- dimensional model space. From these arguments we 
can see that the overbinding in this very limited model space is related to the fact 
that the matrix elements Gj,t(uj' = Ejt) are very attractive for the rather positive 
values of _E ,i an d E lt0 . 

Note that in the limit of HQ = 0, the one- dimensional model-space result for the 
g.s. energy of 4 He is 

E gs = 3(0 - 2.2246) = -6.6738 MeV (for hQ = 0), (9) 

where -2.2246 MeV and are the lowest eigenenergies of the two-body system with 
J=l, T=0 (deuteron) and J=0, T=l, respectively. This (-6.6738 MeV) is the limit 
that the U N=0" curve in Fig. 2 will approach as hQ — > 0. 

In the 2hfl model space, the results for the binding energy of 4 He are reduced 
considerably (see Fig. 3) although they are still larger than the experimental value as 
well as the more exact theoretical value for the RSC potential. Note that the value of 
HQ at which the lowest g.s. energy is obtained is between 16 MeV and 18 MeV in the 
2hQ model space (N=2). This is quite different from corresponding value of hQ=22 
MeV with the 0MI space (JV=0). At htt=l6 MeV, the 2hVt result for the binding 
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energy is 33.6 MeV, which still overbinds the g.s. by a large amount. The reduction of 
the calculated energy with increasing model space can easily be understood from the 
following observations: If the model space is increased, the Pauli operator Q in the 
Bethe-Goldstone Eq.([TJ) ensures that the energy uj\ of the lowest pole in G is shifted 
to higher energies, as this energy correspond to the energy of the lowest 2 particle 
state outside the model space. Therefore the matrix elements of G calculated at the 
same starting energy are less attractive for the larger model space. The net effect of 
enlarging the model space with the appropriately recalculated effective interaction at 
the two-particle level is to reduce the overbinding. 

The results continue to improve as we increase the model space from 2hQ to 4MI. 
Now the lowest g.s. energy of about -31.9 MeV is found when HQ falls between 14 and 
16 MeV. This lowest energy is about 1.7 and 11.4 MeV higher than the corresponding 
values of -33.6 MeV and -43.3 MeV for the 2hfl and OhQ spaces, respectively. 

The dependence of the results on the %Vt value also weakens substantially as we 
go from 2hQ to 2hQ and to 4hQ. We can quantify this dependence by defining a 
dimensionless parameter 

, \E m (Ml + 2) - E„(Q)\ , s 

C a = Average of — 10 

6 2MeV v ; 

that characterizes the "average curvature". For the N=0 curve in Fig. 3, C a is 0.61. 

It decreases to 0.39 for the N=2 curve; and it further reduces to 0.27 for the N=4 

curve. Ultimately, when an infinite Hilbert space is used, the results for the g.s. energy 

should show a complete independence of HQ (i.e., C a defined above vanishes) and 



should converge to the exact result of —24.55 MeV ||14|| . This is indeed the trend we 
are seeing in the OhQ, 2hQ and 4hQ calculations but there is still a considerable gap 
between the 4hfl results and the converged, exact value. Crude extrapolation of the 
C a values for iV=0, 2 and 4 indicates that one may have to do an 8hfl calculation in 
order to reduce C a to less than 0.1, at which point a 10 MeV change in hfl will, on 
average, results in less than 1 MeV change in the ground-state energy. 
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4 Auxiliary single-particle potential 



The discrepancy between the energies obtained in the shell-model calculations of the 
preceding section and the exact result obtained for the RSC potential is due to the 
fact that some effective three- and four-body forces are ignored in our calculations. 

Since shell-model calculations with inclusion of many-body forces are rather in- 
volved ||, we would prefer to find a way to diminish the effect of these many-body 
terms. For that purpose we consider the lowest-order contribution to the three-body 
folded diagrams displayed in Fig. lb. The contribution of this effective three-body 
force to the binding energy within the OhQ model space is represented by the dia- 
gram displayed in Fig.lc. If we introduce an auxiliary s.p. potential for all active 
states of the model space according to the BHF choice 

efHF = u + £ (2 ^ ( ri) +i) <( ' oa)|gr(h; =ti+ e °* )|( * os)) ^ ' (n) 

the contribution of the three-body term of Fig.lc and many higher-order diagrams 
originating from folding would be canceled by corresponding diagrams with s.p. inser- 
tions. Therefore, in this section, we would like to discuss the influence of introducing 
this auxiliary potential on the binding energy calculated in a OhQ model space. 

It should be noted that Eq.flnj) defines an auxiliary potential for the states within 
the model space. For s.p. states outside the model space we assume pure kinetic 
energy. This means that the Bethe-Goldstone Eq.([T]) has been solved with the kinetic 
energy (hi = U) for the s.p. spectrum of the s.p. states outside the model space. 
The resulting G matrix is denoted by G T , as we already did in Eq.(|ll]). When one 
comes to calculate the BHF s.p. energies as defined in Eq.flTTD, this new choice for the 
intermediate energy spectrum is preferred, since our previous choice involves an un- 
specified shift Vq in the one-body potential, which makes it difficult to unambiguously 
relate the starting energy with the BHF s.p. energies in Eq. (|TT|) . 

In order to appreciate the significant effects of using the auxiliary potential, which 
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we will soon discuss, we need to separate the influence of the new choice for the 
intermediate spectrum on the results. To this end, we first evaluate the binding 
energy of 4 He in the OhQ model space without assuming an auxiliary potential. This 
means that we first solve Eq.(|7|) with hi = t{ for the new G- matrix 

h +t 2 + G t {uj = E JjT )] (j)j >T = E, hT <j)j :T , (12) 

so as to determine the matrix elements 

(0s 2 \v^\0s 2 } JjT = (0s 2 \G t (uj = Ej >T )\0s 2 ) J>T . (13) 

We then evaluate the energy according to Eq.(^) and add the Coulomb repulsion 
between the two protons. Results for this calculation without an auxiliary potential 
are displayed in Fig.4 (solid line - "Without aux. pot."). These results are essentially 
the same as those we obtained in the previous section (Fig. 3), indicating that the 
results are rather insensitive to the choice of the intermediate energy spectrum in 



the Bethe-Goldstone equation. This is also consistent with an earlier study [16] in 
which the role of the single-particle potential was examined in some detail at the 
level of two-body effective-interaction calculations. Here we extend those results 
to the case of starting-energy-independent effective two-body interactions. This is 
important since, by eliminating the issue of whether the single-particle insertions 
on intermediate particle lines are responsible for overbinding, we are then forced to 
consider the effective many-body forces that are addressed in the present study. 

We now consider the BHF choice for the s.p. potential. We again use the one- 
dimensional model space for simplicity. In analogy to Eq.(^), we now have 

2e* HF + (0s 2 \G t (lj = E JiT )\0s 2 )j, T ] JjT = E J)T <f> J>T , (14) 

which determines the matrix elements of according to Eq . (|i~3|) . Results for the 



binding energy of 4 He with this choice of the auxiliary potential are also displayed in 
Fig.4 by the solid line labeled 'With BHF pot.". One observes that the use of the 
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BHF auxiliary potential, which has been introduced to measure the effects of terms 
like the three-body folded diagrams of Fig. lb, reduces the calculated binding energy 
drastically Even the minimal value for the binding energy (E gs =-17 .9 MeV), which 
is obtained around HQ = 16 MeV, is well above the "exact result" of Ref. ||14j|. 

For the sake of comparison, we also show in Fig. 4 the results (dotted line) for the 
BHF approximation (restricting the s.p. wave functions to the HO wave functions for 
a given hfl), which can be obtained by replacing Ejt — > 2cq s hf . This approximation 
would correspond to a Otku no-core calculation, assuming the BHF auxiliary potential 
but ignoring the effects of two-body folded diagrams. We see that the two-body 
folded diagrams yield a repulsion of about 2 MeV. Such a repulsive effect has also 
been observed in shell-model calculations within the sd shell [|J . 

5 Conclusions 

(2) 

Energy-independent effective two-body interactions are determined to calculate 
the g.s. energy of 4 He in three no-core HO model spaces with hQ ranging from 10 to 
28 MeV. The results overbind the g.s. of 4 He for all the three spaces and for all the 
hfl values that we have used. The amount of overbinding is largest in the Ohfl space 
and decreases as we increase the size of the model space from OhQ to 2hQ and to 4M7. 
The dependence of the calculated g.s. energy on the value of HQ also becomes weaker 
as the size of the model space increases. However, even in the 4hQ calculation, the 

(2) 

lowest g.s. energy of -31.9 MeV, which we obtained with vl s at hQ=16 MeV, still 
overbinds 4 He by 3.6 MeV when compared with the experimental binding energy and 
by 7.3 MeV when compared with the binding energy obtained in the GFMC approach 

0- 

This overbinding is caused by the fact that we are neglecting some effective many- 
body forces. Our results using a BHF auxiliary potential, which includes certain 
effective many-body-force terms, show significant effects on no-core shell-model cal- 
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culations, especially when the model space is small. More studies are required to find 
an optimal auxiliary potential, which minimizes the effects of such many-body forces. 
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Figure captions 



Fig.l Contributions to the effective three-body force originating from the so-called 
Q-box (a) and folded diagrams (b). The diagram (c) represents the contribution of 
diagram(b) to the ground-state energy within a Ok£l model space. The bare lines refer 
to s.p. states within the model space while the "railed" lines respresent s.p. states 
outside the model space. 

Fig. 2 The two-body matrix elements, (0s 2 \G(uj')\0s 2 }j yT (solid lines) and (0s 2 \v^ |0s 2 ) JyT 
(dotted lines) for the 4MI model space with TiVt=lQ MeV as a function of the starting 

(2) 

energy u'. Note that the TBMEs of are independent of cu'. 

Fig. 3 The ground-state energy of 4 He obtained from function of TiVt and 

N, the size of the model space. 

Fig. 4 The g.s. energy of 4 He as a function of hfl obtained from the Oh£l calculation. 
Results are represented for using kinetic energies (labeled: Without aux. pot.) or 
BHF s.p. energies (label: With BHF pot.) as an auxiliary potential. For a comparison 
the dashed line displays the energies obtained with the BHF approximation. 
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